///- 


IMPACT OF TROPOSPHERIC AEROSOL ABSORPTION 
ON OZONE RETRIEVAL FROM BUY MEASUREMENTS 


O. Torres' and P.K. Bhartia 2 
Raytheon STX Corporation 
4400 Forbes Blvd, Lanham, Md, 20706 


(Submitted to Journal of Geophysical Research, Atmospheres) 


October 27, 1998 


'Current affiliation: Joint Center for Earth Science Technology, University of MD, 
Baltimore, MD 21250 


2 NASA/Goddard Space Flight Center, Code 916, Greenbelt, MD, 20771 



Abstract 


The impact of tropospheric aerosols on the retrieval of column ozone amounts using 
space borne measurements of backscattered ultraviolet radiation is examined. Using radiative 
transfer calculations, we show that uv-absorbing desert dust may introduce errors as large as 
10% in ozone column amount, depending on the aerosol layer height and optical depth. Smaller 
errors are produced by carbonaceous aerosols that result from biomass burning. Though the error 
is produced by complex interactions between ozone absorption (both stratospheric and 
tropospheric), aerosol scattering, and aerosol absorption, a surprisingly simple correction 
procedure reduces the error to about 1%, for a variety of aerosols and for a wide range of aerosol 
loading. Comparison of the corrected TOMS data with operational data indicates that though the 
zonal mean total ozone derived from TOMS are not significantly affected by these errors, 
localized affects in the tropics can be large enough to seriously affect the studies of tropospheric 
ozone that are currently undergoing using the TOMS data. 

Introduction 

The inversion of satellite measured backscattered ultraviolet (buv) radiances has proven 
to be an extremely reliable method of globally mapping the spatial and temporal distribution of 
the atmospheric ozone content. Since the first space borne buv sensor was launched in 1970, 
several ozone sensing experiments (TOMS, buv, GAME) have used observations of the buv 
radiation to continuously monitor the stratospheric ozone layer. Most future satellite missions 
(NPOESS-OMPS, EOS- CHEM, ENVISAT, METOP ) to study the ozone layer are expected to 
include buv-based ozone sensors. It was recognized very early [Dave, 1978] that the scattering 
and absorption effects of aerosols (both tropospheric and stratospheric) modify the ultraviolet 



radiation field and affect the retrieval of total ozone content and its vertical distribution. Though 
the impact of background stratospheric aerosols on the ozone retrieval is small, large volcanic 
eruptions such as El Chichon (1982) and Mt. Pinatubo (1991) produce easily recognizable 
artifacts in the total ozone field [Bhartia et. al, 1993, Torres et. al., 1995] and in the retrieved 
vertical distribution [Mergenthaler, 1985; Torres and Bhartia, 1995]. 

This paper discusses the effect that tropospheric aerosols have on the TOMS (Total 
Ozone Mapping Spectrometer) ozone retrieval process. Though earlier studies [e.g., Dave, 1978] 
looked at this problem, they used a limited number of synthetic aerosol models that were not tied 
to particular types of aerosols found in the Earth's atmosphere. A fresh look at the problem has 
become important for two reasons: 1) a new method of mapping aerosols using TOMS 
measurements has become available [Torres et al, 1998; Herman et al., 1997], and 2) residual 
techniques applied to the TOMS data are now also utilized to infer tropospheric ozone content 
[. Fishman et al., 1990; Hudson and Thompson, 1998; Ziemke et al., 1998]. Since tropospheric 
ozone is about 10% of the total ozone, a small error in total ozone can propagate to the derived 
tropospheric ozone yielding a much larger percent error in the derived tropospheric ozone 
content. Moreover, most of the studies of tropospheric ozone focus on precisely the areas where 
the atmosphere has a large aerosol loading due to biomass burning or other pollution episodes. 

We begin with an analysis of the effect of tropospheric aerosol types on the buv 
radiances in section 2. The resulting errors in retrieved ozone are examined in section 3 by using 
synthetic aerosol contaminated radiances in the TOMS algorithm. Finally a method to correct for 
the tropospheric aerosol effect is discussed. 



2. Tropospheric aerosol effect on backscattered radiances 
-Aerosol Models 

A set of aerosol models representative of the most abundant tropospheric aerosol types 
were selected for this analysis. Table 1 summarizes the optical properties of the aerosol models 
[Torres et. al., 1998], Non-absorbing tropospheric aerosols of anthropogenic origin are 
represented by a sulfate aerosol model [Shaw, 1979] denoted by S in Table 1. Model Cl 
represents weakly UV-absorbing carbonaceous particulate representative of aerosol material from 
the smoldering phase of combustion. It is the most abundant aerosol type in regions where 
intensive biomass burning takes place on a seasonal basis. Model C2 is representative of the 
more absorbing carbonaceous aerosol particles associated with the flaming phase of the fire 
generally found near the source region. Since data on the refractive index of carbonaceous 
aerosols in the UV are not available, its spectral dependance is not taken into account, and the 
assumed values in Table 1 are based on measurements in the visible [Patterson and McMahon , 
1984]. Model D represents strongly absorbing desert dust aerosols. The spectrally dependent 
refractive index of the D model is taken from measurements of Saharan desert dust absorption 
by Patterson [1977]. 


Table 1. Aerosol Models used in the simulation study 


Aerosol Model 

Parameters 

Refractive Index 

S 

0.07 2.03 

1.43 - O.OOOi 

Cl 

0.14 1.45 

1.55 - 0.01 5i 

C2 

0.14 1.45 

to 

m 

o 

o 

1 

to 

to 

D 

0.25 2.20 

1.57 - 0.015i 




The aerosol vertical distribution of the absorbing aerosols is assumed to be a gaussian 
distribution characterized by the height of maximum aerosol concentration (peak of the 
distribution) and its half-width. The non-absorbing aerosol vertical distribution is assumed 
maximum at the surface and decreasing exponentially with height. Calculations were done for 
several values of aerosol optical depth and aerosol-layer altitude. Low latitude ozone profiles 
[McPeters et al, 1996] with integrated ozone contents of 275 and 325 DU were assumed. The 
effect of each aerosol model was analyzed separately i.e. no aerosol mixtures were considered. 

Synthetic aerosol contaminated radiances were used to simulate the intensities measured 
from space when the atmosphere contains suspended particulate represented by the aerosol 
models in Table 1 for an optical depth of 1 .0 (at 380 nm). For the sulfate aerosol model (S 1 ), an 
optical depth of 1 .0 should be regarded as representative of the most severe cases of tropospheric 
sulfate contamination. For carbonaceous and mineral aerosol models, this value of optical depth 
can be considered typical. Optical depths of 1 .0 and larger are commonly observed during the 
large scale biomass burning events known to take place seasonally in South America [Gleason et. 
al., 1998] and Equatorial Africa. Sunphotometer measurements of desert dust aerosols 
periodically report values in excess of 1 .0 during dust outbreaks [ Pinker et al . , 1 994; 

D'Almeida, 1987; Faizoun et al., 1994]. 

- Spectral dependance of the aerosol effect 

In addition to aerosol physical properties (particle size distribution and complex refractive 
index), the change in upwelling ultraviolet radiance at the top of an aerosol laden atmosphere 



also depends on the aerosol optical depth, its location in the atmosphere and reflectivity of the 
underlying surface. The change in radiance, with respect to a molecular atmosphere with a 
vertically distributed ozone amount, is shown in Figure 1 for the aerosol models used in this 
work (solid lines). The transition from a radiance increase at all wavelengths for the 
non-absorbing sulfate-aerosol model to a radiance decrease for the highly absorbing mineral-dust 
model is produced by the competing effects of scattering and absorption by the particles. A 
detailed discussion of the dependance of the aerosol effect on other factors such as altitude of the 
aerosol layer and reflectivity of the underlying surface is presented by Torres et al. [1998]. 

The aerosol effect in the 340-380 nm range is approximately linear with wavelength for 
all aerosol models. The slight non-linearity of the D model is related to the spectrally dependent 
imaginary component of the refractive index used in the calculations. A clear departure from 
linearity is observed at wavelengths shorter than about 340 nm. This deviation from linearity is 
not directly associated with the aerosol spectral optical properties. Rather, it is the result of the 
complex interaction of the scattering and absorption (by ozone and aerosols) processes of both 
Rayleigh and Mie scattered photons. This is clearly illustrated by the dashed lines in Figure 1 , 
which show that when the atmosphere does not contain ozone, the linearity of the aerosol effect 
extends to the 310-340 spectral range. 

3. Aerosol Impact on TOMS ozone retrieval 

In this section we present quantitative estimates of the errors introduced by aerosols in the 
ozone amounts retrieved from satellite measurements of buv radiation using the TOMS 
algorithm. A very short and simple description of the TOMS algorithm follows. For a more 



detailed discussion of the TOMS ozone retrieval process, the reader is advised to see the TOMS 
User's Guide [ McPeters et al . , 1 996] . 


- TOMS Algorithm 

The current operational (version 7) TOMS total-ozone algorithm is essentially a 2-step 
process. In the first step, an initial estimate of total ozone is obtained by assuming a "Rayleigh 
model" of the atmosphere. In this model one approximates the true atmosphere by a a molecular 
atmosphere containing a vertically distributed ozone amount bounded at the bottom by a 
Lambertian reflecting surface. For clear-sky conditions, or for fully-cloudy scenes containing 
optically thick clouds, the reflectivity of the surface is adjusted to match the measured radiances 
at one wavelength (380 or 360 nm), called the reflectivity wavelength (A, R ). For scenes containing 
broken clouds, radiances are obtained by linearly mixing the cloudy and clear scenes (using 
nominal values of surface and cloud reflectivities) to match the radiances at A. R . (For further 
details the readers are referred to the relevant TOMS User guides.) The algorithm predicts the 
radiances at other TOMS wavelengths by detailed radiative transfer computations that account 
for ozone absorption, multiple scattering and polarization effects, Earth’s sphericity effects, 
surface reflection, as well as effects of instrument bandpass and atmospheric temperature. A 
simple correction is applied for Rotational Raman Scattering (Ring effect) using the model 
developed by Joiner et al. [1995]. By interpolation of pre-computed tables that vary with latitude 
and total ozone amount, the algorithm finds a value of total ozone that explains the ratio of 
measured radiances at a pair of wavelengths (3 17.5/33 1 .2 nm). 



Careful analysis of TOMS data at weakly ozone absorbing wavelengths (>330 nm) 
indicated that though the Rayleigh model (RM) worked surprisingly well in most cases, there 
were some notable exceptions. These included highly non-lambertian surfaces (e.g., sea glint) 
and UV-absorbing aerosols (volcanic ash, smoke, and mineral dust). Careful study of these 
deviations from the RM has led to the discovery of a fundamental new technique for detecting 
tropospheric aerosol from space [Herman et al., 1997, Torres et al., 1998], To account for these 
effects (as well as for the selection of a better ozone profile, which will not be discussed here), 
the Version-7 TOMS algorithm includes a second step. In this step the difference in the 
logarithm of the measured radiances and those calculated from the RM, called a residue, is 
assumed to be linear with wavelength. (Note that, by definition of the RM, the A R residue is 
zero). Ozone is then recalculated by applying this correction, using a first order Taylor series 
expansion around the previously calculated ozone value. (For detail see the references cited 
above.) The key point to note is that for this correction to be accurate, the residuals must be 
linear with wavelength. We examine this assumption next. 

- Simulation experiment 

A simulation study was conducted in order to quantify the error in the retrieved ozone- 
column amount when the atmosphere is loaded with aerosol particles. Radiative-transfer 
calculations at the TOMS wavelengths were performed to generate the aerosol contaminated 
radiances emanating at the top of an atmosphere of known total ozone content. Calculations 
were carried out for the aerosol models discussed in section 2. To evaluate the aerosol effect 
under different conditions, radiative transfer calculations were done for several values of the 380 
nm aerosol optical depth (0.5, 1., 2. and 4.), and different locations of the absorbing aerosol layer 



(0.5, 1.5, 3. and 6 km) above the ground. The aerosol and ozone vertical distribution were 
discussed in section 2. 


The aerosol contaminated radiances were input to the TOMS algorithm. The difference 
between the retrieved ozone amount and the actual ozone content is the retrieval error associated 
with the aerosol interference. The error is expressed in percent of the total ozone amount. 
Calculations using different total ozone amounts indicate that although the absolute error 
increases with increasing total ozone amount, the relative or percent error does not depend on the 
total-ozone content provided that the tropospheric component does not change. The implication 
of tropospheric ozone variability will be discussed in a subsequent section. 

Since, in the simulation analysis the true atmospheric-ozone content is known, the 
spectral-dependance departure (from a molecular atmosphere) shown in Figure 2 (dashed lines) 
in terms of residues is computed. The results reflect the interaction of the several radiative- 
transfer processes taking place in the atmosphere: molecular and particle scattering as well as 
absorption by ozone molecules and aerosol particles of both Rayleigh- and Mie-scattered 
radiation. Also shown, are the residues predicted by the TOMS algorithm at 312.5, 331.2 and 
380 nm. The straight-solid lines are obtained by assuming a linear dependance of the residues 
in the 331-380 range and extrapolation to 312.5 nm. Although the linear function is a good 
representation of the wavelength dependance of the residues produced by weakly absorbing 
aerosols (model Cl), there is a significant departure from linearity for highly absorbing and 
non-absorbing aerosols. Because of this, when the atmosphere is loaded with highly absorbing 
aerosols, the retrieved ozone amount is in error. 



- Angular dependance of the aerosol effect 

The simulated ozone-retrieval error for the aerosol models used in this analysis are 
shown in Figure 3 for typical TOMS viewing conditions. In general, the effect of the non- 
absorbing aerosol particles (model SI) is to induce a scan-angle dependent ozone 
over-estimation of up to about 2%. As the satellite-zenith angle increases beyond 50%, the 
error decreases, becoming negative (-2%) at the extreme TOMS scan angles. The angular 
dependance is associated with the aerosol-scattering phase function. The angular dependance of 
the carbonaceous-aerosol models (Cl and C2) is not as pronounced, whereas the mineral aerosol 
model (D) shows a well defined scan angle variability. 

- Effect of aerosol absorption 

The error in the retrieved ozone amount in the presence of absorbing particles increases 
with increasing aerosol optical depth and height of the absorbing aerosol layer. Contours 
illustrating the resulting ozone artifact as a function of aerosol optical depth and layer height are 
shown in Figure 4. The sign of the ozone error resulting from the presence of weakly-absorbing 
particles (model Cl) depends on the combined effect of the aerosol-layer optical depth and its 
altitude above the surface. As shown in Figure 4a, small (less than 1%) positive errors are 
produced by aerosol layers located as high as 2 km above the surface and with optical depths as 
large as 2.0 . Aerosol layers higher than about 3 km, or optical depths larger than about 3, are 
required to produce negative errors. Thus, for realistic conditions of aerosol optical depth as 
large as 3 and aerosol layer location of 3 km, the ozone error due to Cl aerosols is about 1%. 



The more absorbing carbonaceous aerosol model (C2) produces predominantly negative 
errors except when the aerosol layer lies a few hundred meters above the surface, as shown in 
Figure 4b. Ozone underestimations between 3% and 4% are possible for aerosol layers at 3 km 
and optical depth 3 or larger. The effect of strongly absorbing mineral aerosols (model D) on the 
ozone retrieval is shown by the contour plot in Figure 4c. This aerosol type yields negative 
ozone-amount errors regardless of aerosol optical depth and location in the atmosphere. For 
typical optical depth values between 1 and 2 and aerosol-layer height between 3 and 5 km, 
ozone-amount errors as large as 10% are possible. 

- Error sources 

In this section we examine the causes of the aerosol errors we have just discussed. From 
Figure 1 we already know that the linear model assumed by the TOMS-V7 algorithm is correct in 
absence of atmospheric ozone, but it becomes inaccurate when ozone is added. To examine if the 
effect is primarily due to change in effective tropospheric airmass factor (photon path length) due 
to interaction between tropospheric ozone and aerosols, we ran our mineral dust model with and 
without tropospheric ozone. These results, shown in Figure 5, indicate that 50-75% of the ozone 
error remains even in the absence of tropospheric ozone, even for aerosols close to the surface . 
Since aerosols so low in the atmosphere cannot directly alter the path of photons in the 
stratosphere, the error results from the fact that ozone reduces the number of photons reaching 
the troposphere and hence interacting with aerosols. The remaining error (shown by dotted line), 
attributable to tropospheric ozone, is probably due to change in the path of the photons in the 
troposphere due to aerosols. Part of this simply is a shielding effect, i.e., a highly absorbing layer 
would partially shield the ozone below it. However, the fact that the dotted line doesn't go to 



zero for very low altitude aerosols, indicates that the absorption of photons by ozone above the 
aerosol layer is also reduced. This is probably because there is less multiple scattering in 
presence absorbing aerosols than in the Rayleigh atmosphere. 

The foregoing discussion highlights the fact that the effect of aerosols on buv radiation is 
strongly affected by atmospheric ozone, and hence cannot be corrected by any form of monotonic 
correction scheme, linear or otherwise. This has implication for other methods of retrieving 
ozone from buv radiances, e.g., the DOAS method used in processing the GAME buv data 
[European Space Agency, 1995]. 

4. Correction Procedure 

Because the aerosol effect on the buv radiation depends on several different parameters 
(complex refractive index, particle size distribution, aerosol-optical depth and aerosol-layer 
height), it is extremely difficult to accurately correct the measured buv radiances for the radiative 
transfer effects of absorbing aerosols. The difficulty lies in the fact that required information on 
the aerosol properties and its vertical distribution is generally unavailable. 

To search for a simple correction scheme, we examined the relationship between the 
ozone error and the Aerosol Index (AI), which is simply the 340 nm (33 1 nm for EP/TOMS) 
residue in N-value units. (If the residues are defined using the natural logarithm, N-value, 
obtained by multiplying it with -100 log I0 (e), converts it into a more convenient unit). Based on 
previous work [Herman et ai, 1997] we know that AI provides a convenient way of tracking 
absorbing aerosols all over the globe. Figure 6 shows the ozone error for the aerosol models 



used in this analysis as a function of AI. As illustrated in Figure 6a for aerosol model Cl, at a 
fixed aerosol-layer height, both the ozone-retrieval error and the magnitude of the AI increase 
with increasing aerosol optical depth are shown. The slope of the ozone error - AI relationship 
shows only a slight variation for different heights of the aerosol layer. Results of similar analysis 
for aerosol models C2 and D are illustrated in Figures 6b and 6c respectively. In spite of the 
somewhat larger variability of the ozone error - AI relationship with aerosol layer height for 
model D, the observed relationship is very similar for all the aerosol models. 

For other viewing geometries, the AJ-error relationship for the absorbing aerosol models 
does not deviate significantly from the linear function in Figure 6. For the non-absorbing 
aerosol model , however, the Al-error relationship does vary significantly with geometry. As 
shown in Figure 3, the ozone error can be either negative or positive depending on the geometry 
of the observation. The variability associated with aerosol type, optical depth, aerosol layer 
height and viewing geometry is shown on the scatter plot in Figure 7. From this analysis, it is 
apparent that the relationship between the retrieval error and AI for the absorbing aerosol models 
is general enough over a wide range of conditions to be described mathematically by a simple 
linear function as shown by the linear fit to the data in Figure 7. 

The well-defined relationship between the positive values of the AI and the aerosol- 
absorption artifact in the retrieved ozone amount, can be used to apply a first-order correction to 
the ozone product. This correction method does not require any information on the aerosol 
microphysical properties or location in the atmosphere. An additional advantage is that such a 
correction can easily be applied to the retrieved ozone amounts without requiring knowledge of 



the specific viewing geometry. 


Based on the previous discussion, the ozone artifact (in Dobson units) associated with the 
presence of absorbing aerosols is 

£ = O.OlkAIQ ret (1) 

where Q rct is the retrieved ozone amount, AI is the positive aerosol index, and k is the slope of 
the linear fit in Figure 9. The value of k is 1 . 12 for the 340-380 AI and 1 .20 for the 331-360 AI 
(used for EP/TOMS). 

5. Aerosol Effect on Nimbus7-TOMS data 

Validation of the above described relationship in the TOMS data is not direct, since the 
ozone variability associated with the aerosol interference must be separated from the natural 
spatial and temporal variability of the ozone field. The natural variability may include 
tropospheric ozone increases as a biomass burning by-product, or possible reductions as result of 
ozone processing on the surface of mineral dust aerosols [ Dentener et. al, 1996]. 

A qualitative verification of the relationship between the absorbing aerosol induced ozone 
underestimation and the AI is observed over the Saharan Desert. The effect is more easily 
observed during the Summer when the aerosol layer is generally higher (about 4km or more 
above the surface) yielding larger AI values for the same amount of aerosol. These regions of 
high AI are frequently associated with total ozone amounts lower than in neighboring areas 
where the aerosol presence, as detected by the AI, is not as strong. 



Over the biomass burning areas, on the other hand, the verification of the theoretically 
derived AI - ozone-error relationship is more complicated. Although some areas of reduced 
ozone and large AI can be found, there are areas where the ozone content is enhanced in relation 
to the surrounding areas in spite of large AI values. A likely explanation of this reverse 
correlation is that the natural ozone increase generally associated with biomass burning events 
could be offsetting the aerosol absorption artifact. 

Plate 1 shows the global geographical distribution of the error in the TOMS derived 
ozone field on August 23, 1985, as obtained using equation 1. As expected, the retrieval errors 
are localized over the arid regions of the world and over areas known to be affected by large scale 
biomass burning such as South America and Central Africa. On this particular day, derived errors 
as large as 5% (about 14 DU) are observed over the western Saharan Desert, the coast of Angola, 
and Brazil. The effect of the aerosol artifact on the zonally averaged total ozone data is 
negligible (less than 1%). 

6. Summary and Conclusions 

The radiative transfer processes of aerosol scattering and absorption affect the retrieval of 
ozone-column amount by the TOMS sensor. The algorithm-predicted linear spectral dependance 
of the effect of non-absorbing and weakly absorbing particles on the backscattered radiances is 
very close to the observed wavelength dependance. For this reason, the ozone retrieval error in 
the presence of non-absorbing sulfate and weakly-absorbing carbonaceous aerosols is small. The 
deviation of the observed spectral dependance from the algorithmic assumption is larger when the 
atmosphere is loaded with moderately absorbing carbonaceous particulate and highly-absorbing 



mineral dust. This causes significantly larger errors to occur. 


The absorbing-aerosol induced error consists of a pseudo-reduction of the atmospheric 
ozone column over those areas where large amounts of absorbing aerosols are present at 2 km or 
higher in the atmosphere. The retrieval error increases with increasing aerosol optical depth and 
aerosol-layer height. The ozone under-estimation due to the effect of biomass burning aerosols are 
generally small (2 % or less), since smoke layers generally reside in the lowest 2 km of the 
atmosphere. Only when the smoke layer rises significantly or is present in large amounts, are 
errors larger than 2% expected. On the other hand, errors as large as 8% are produced by layers 
of desert dust at 4km or higher. 

The close linear relationship between the ozone error associated with absorbing aerosols 
and the TOMS Aerosol Index is the most important finding of this work. It has been shown that 
regardless of the aerosol type, optical depth, aerosol layer height and viewing geometry, there 
exists a well defined linear relationship between the error (as a percent of the total column 
amount) and the value of the Aerosol Index. This relationship is general enough to allow the 
application of a first-order correction to the ozone data. Because the variability associated with 
the viewing geometry is small, the correction can be applied directly to the level-3 griddled data. 
The corrected ozone data is within 1 % of the true value. 

Correcting the ozone retrieval for the aerosol effect reduces the uncertainty of the 
tropospheric ozone content derived as a difference between the TOMS total ozone and other 
independent estimates of the stratospheric ozone content. The use of uncorrected TOMS total- 



ozone amounts, in the application of the residual techniques, results in spurious tropospheric 
ozone deficits. 
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Figure Captions 

Figure 1 . Spectral dependance of the change in backscattered radiance, with respect to a 

molecular atmosphere with ozone absorption for the four aerosol models in Table 1 (solid 
lines). The obtained spectral dependance when the atmosphere does not contain ozone is 
shown as the dashed lines. The symbols indicate the wavelengths of the calculations. The 
aerosol optical depth (at 380 nm) is set to unity and the absorbing aerosol layer is placed 
at 3.0 km. Other assumed conditions are: ozone content , 275 DU; reflectivity of the 
surface, 0.05; solar zenith angle, 20%, nadir viewing. 

Figure 2. Residues (see text) as function of wavelength for the aerosol models used in the 

analysis (dotted lines). The solid straight lines represent the algorithm assumed spectral 
dependance of the residual quantities. All modeling conditions are the same as in figure 1 . 

Figure 3. Percent error in retrieved total ozone as a function of satellite zenith angle, for the set 
of assumed aerosol models at relative azimuth angles typical of a sun-synchronous orbit at 
low latitudes. Other modeling conditions as in figure 1. 



Figure 4. Contours describing the ozone error (in percent) introduced by absorbing aerosols as a 
function of the aerosol optical depth and aerosol layer height above the ground. Different 
panels correspond to: a) weakly absorbing carbonaceous aerosols (model Cl); b) 
Moderately absorbing carbonaceous aerosol (model C2); c) Mineral dust (model D). 

Other modeling conditions as in figure 1 . 

Figure 5. Percent ozone retrieval error as a function of aerosol layer height for the mineral aerosol 
model (D) of optical depth 1.0 (dashed line). The dotted line represents the component of 
the total error associated with the 'ozone-masking' effect of the aerosol layer (see text). 

The thin solid line represents the part of the total error resulting from the error in 
predicting the spectral dependance of the aerosol effect on the backscattered radiances. 

Figure 6. Relationship between the percent ozone error and the TOMS measured Aerosol Index 
for (a) the Cl model, (b) the C2 model and (c) the D model. The lines represent different 
heights of the aerosol layer: solid line, 0.5 km; dotted line, 1.5 km; dashed line, 3 km and 
dot-dashed line, 6 km. The symbols correspond to calculations for optical depth values of 
0, 0.5,1, 2 and 4, increasing in the direction of the arrows. 

Figure 7. Scatter plot illustrating the variability of the ozone error-Aerosol index relationship for 
simultaneously changing aerosol type (C 1 , C2, D); aerosol layer height (0.5, 1.5,3 and 6 
km); aerosol optical depth (0.5, 1.0, 2.0 and 4.0); relative azimuth angle (0% to 180%, 
30% steps) and satellite zenith angle (0% to 64%, 2% steps). Also shown is the linear fit to 


the synthetic data points. 



Plate 1. Global Geographical Distribution of the error in Nimbus7-TOMS retrieved total ozone 
amount on August 23, 1985, as derived by making use of the TOMS measured Aerosol 
Index. 
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